Observation of Ultra-high-energy Cosmic Rays 
with the ANITA Balloon-borne Radio Interferometer 



O 
(N 



(N 



6 

(N 

> 

in 

en 

o 
p 

o 
o 



% 



S. Hoover', J. Nam^, P. W. Gol-ham^ E. Grashorn'', P. Allison^ S. W. Barwick^ J. J. Beatty'', K. Belov', 

D. Z. Besson*', W. R. Binns^, C. Chen^ P Chen^'', J. M. Clem^, A. Connolly'", P F. Dowkontt'', 
M. A. DuVemois^'2, R. C. Field", D. Goldstein^ A. G. Vieregg', C. Hast", C. L. Hebert\ M. H. Israel^, 

A. Javaid^. J. KowalskP, J. G. Learned^ K. M. Liewer'^ J. T. Link^''*, E. Lusczek'^, S. Matsuno^ 
B. C. Mercurio", C. Miki\ P Miocinovic^ C. J. Naudet'^ J. Ng", R. J. Nichol'°, K. Palladino'^, K. Reil", 
A. Romero-Wolf^ M. Rosen^ L. Ruckman^ D. Saltzberg', D. Seckel*^, G. S. Varner^ D. Walz", F Wu^' 

Dept. of Physics and Astronomy, Univ. of California, Los Angeles, 

CA 90095. ^Dept. of Physics, Ewha Womans University, Seoul, 

South Korea. Dept. of Physics and Astronomy, Univ. of Hawaii, 

Manoa, HI 96822. ^Dept. of Physics, Ohio State Univ., Columbus, 

OH 43210. ^Dept. of Physics, Univ. of California, Irvine, 

CA 92697. ^Dept. of Physics and Astronomy, Univ. of Kansas, 

Lawrence, KS 66045. Dept. of Physics, Washington Univ. in St. Louis, 

MO 63130. ^Dept. of Physics, Grad. Inst, of Astrophys., 

& Leung Center for Cosmology and Particle Astrophysics, National Taiwan University, 

Taipei, Taiwan. Dept. of Physics, Univ. of Delaware, Newark, 

DE 19716. ^^^Dept. of Physics and Astronomy, University College London, 

London, United Kingdom. ''^^SLAC National Accelerator Laboratory, Menlo Park, 

CA, 94025. School of Physics and Astronomy, Univ. of Minnesota, 

Minneapolis, MN 55455. ^^Jet Propulsion Laboratory, Pasadena, 

CA 91109. ^'^Currently at NASA Goddard Space Flight Center, Greenbelt, MD, 20771. 

We report the observation of sixteen cosmic ray events of mean energy of 1.5 x lO'^ eV, via radio pulses 
originating from the interaction of the cosmic ray air shower with the Antarctic geomagnetic field, a process 
known as geosynchrotron emission. We present the first ultra-wideband, far-field measurements of the radio 
spectral density of geosynchrotron emission in the range from 300-1000 MHz. The emission is 100% linearly 
polarized in the plane perpendicular to the projected geomagnetic field. Fourteen of our observed events are 
seen to have a phase-inversion due to reflection of the radio beam off the ice surface, and two additional events 
are seen directly from above the horizon. 



The origin of ultra-high energy cosmic rays (UHECR) re- 
mains a mystery decades after their discovery yj, |2[]. Key 
to the solution will be increased statistics on events of high 
enough energy (> 3 x lO'^ eV) to elucidate the endpoint of 
the UHECR energy spectrum as seen at Earth. The primary 
difficulty is the extreme rarity of events at these energies. 
Despite steady progress with experiments such as the Pierre 
Auger Observatory, there remains room for new methodolo- 
gies. Cosmic rays have been detected for decades via im- 
pulsive radio geosynchrotron emission yl I5l 4l6ll but until 
now not in this crucial energy range, which offers the pos- 
sibility of pointing the UHECRs back to their sources. We 
present data from the Antarctic Impulsive Transient Antenna 
(ANITA) 1.2 1.1 which represents the first entry of radio tech- 
niques into this energy range. We find 16 UHECR events, 
at least 40% of which are above lO'^ eV, and we show com- 
pelling evidence of their origin as geosynchrotron emission 
from cosmic -ray showers. Our results indicate degree-scale 
precision for reconstruction of the UHECR arrival direction, 
lending strong credence to efforts to develop radio geosyn- 
chrotron detection as a competitive method of UHECR parti- 
cle astronomy. 

Geosynchrotron emission arises when the electron-positron 



particle cascade initiated by a primary cosmic ray encounters 
the Lorentz force in the geomagnetic field. The resulting ac- 
celeration deflects the electrons and positrons and th ey begi n 
to spiral in opposite directions around the field lines BlTlllSll . 



In air, the particles' radiation length is of order 40 g cm^ , a 
kilometer or less at the altitudes of air shower maximum de- 
velopment. Particle trajectories form partial arcs around the 
field lines before they lose enough energy to drop out of the 
shower. The meter-scale longitudinal thickness of the shower 
particle 'pancake' is comparable to radio wavelengths below 
several hundred MHz; thus the ensemble behavior of all of the 
cascade particles yields forward-beamed synchrotron emis- 
sion which is partially or fully coherent in the radio regime. 
Therefore, the resulting radio impulse power grows quadrat- 
ically with primary particle energy, and at the highest ener- 
gies, yields radio pulses that are detectable at large distances. 
Current systems under development for detection of these ra- 
dio impulses are co-located with cosmic-ray particle detectors 
on the ground to aid in cross-calibration il4l4l6ll . They de- 
tect showers with primary energies in the lO'^^'^ eV range 
because of their limited acceptance. No such system has re- 
ported a sample of > lO'^ eV UHECR events. 

The ANITA long-duration balloon payload is launched 
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FIG. 1 : An example of interferometric maps of relative correlated 
intensity for both Hpol (top) and Vpol (bottom) from event 3623566 
which occurred in a region of Antarctica where the geomagnetic in- 
clination gave an appreciable Vpol component for the shower radio 
emission. The sidelobes are residuals from the relatively sparse sam- 
pling of the ANITA interferometer baselines. Such maps are used to 
verify the location of the emission source on the Antarctic continent, 
and exclude emission that arises from known anthropogenic sources. 



from Williams Field near McMurdo Station, Antarctica. It 
takes advantage of the stratospheric South Polar Vortex to cir- 
cle the Antarctic continent at altitudes of 35-37 km while syn- 
optically observing an area of ice of order 1 .5M km^. During 
flight, ANITA records all nanosecond-duration radio impulses 
over a 200-1200 MHz radio frequency band. The threshold is 
a few times the received power of thermal emission from the 
ice, ~ 10 picoWatts. The direction of detected signals, deter- 
mined by pulse-phase interferometric mapping (Fig.[Tl[2l|]), 
is localized to an angular ellipse of 0.3° x 0.8° (elevation x 
azimuth) which is projected back onto the continent to deter- 
mine the origin of the pulse. ANITA's mission is the detection 
of ultra-high energy neutrinos via linearly-polarized coherent 
radio Cherenkov pulses from cascades the neutrinos initiate 
within the ice sheets. Virtually all impulsive signals detected 
during a flight are of anthropogenic origin, but such events 
can be rejected with high confidence because of their associa- 
tion with known human activity, which is carefully monitored 
in Antarctica. For its first flight, during the 2006-2007 Aus- 
tral summer, ANITA's trigger system was designed to max- 
imize sensitivity to linearly polarized radio pulses, but pur- 
posely blinded to the plane of polarization. However, the 
entire polarization information - both vertical and horizon- 
tal (Vpol and Hpol) - was recorded for subsequent analysis. 
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FIG. 2: Map of locations of detected UHECR events superimposed 
on a Radarsat image of relative microwave radar backscatter ampli- 
tude of the Antarctic continent. The red diamonds are the reflected- 
event locations, the black squares are the two direct-event locations. 
The dash-dot line indicates the limit of ANITA's field-of-view for the 
flight. Note that the portion of ANITA's field-of-view that includes 
the ocean was always covered by sea-ice during the flight. 



Since radio pulses of neutrino origin strongly favor vertical 
polarization, due to the geometric-optics constraints on the 
radio Cherenkov cone as it refracts through the ice surface, 
we used the Hpol information as a sideband test for our blind 
neutrino analysis. 

Our results were surprising: while the neutrino analysis 
(Vpol) gave a null result, a statisti cally significant sample of 
6 Hpol events was found initially 12211 . and a more sensitive 
analysis now yields 16. These events are randomly distributed 
around ANITA's flight path (Fig.|2]i, uncorrected in location 
to human activity or to each other, but closely correlated to 
each other in their radio pulse profile and frequency spec- 
trum (Fig. [3] Top). Their measured planes of polarization are 
found in every case to be perpendicular to the local geomag- 
netic field (Fig. m, as expected from geosynchrotron radia- 
tion. With two exceptions, the events reconstruct to locations 
on the surface of the ice; the two exceptional cases have di- 
rectional origins above the horizon, but below the horizontal 
(from stratospheric balloon altitudes the horizon is about 6° 
below the horizontal). Earth-orbiting satellites are excluded 
as a possible source since the nanosecond radio temporal co- 
herence observed is impossible to retain for signals that prop- 
agate through the ionospheric plasma, which is highly dis- 
persive in our frequency regime. The fourteen below-horizon 
events are inverted compared to the two above-horizon events, 
as expected for specular reflection (Fig. [3] Top). From these 
observations we conclude that ANITA detects a signal, seen 
in most cases in reflection from the ice sheet surface, which 
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FIG. 3: Top: overlay of the 16 UHECR event Hpol pulse shapes, 
showing the inverted phase for the 14 reflected events (in blue) com- 
pared to the two direct events (in red). Inset: Average pulse pro- 
file for all events. Bottom: Flux density for both the averaged di- 
rect and averaged reflected events. In each case the data are con- 
sistent with an exponential decrease with frequency: the fitted co- 
efficients of decrease with frequency are (180 ± 13 MHz)^ , and 
(197 ± 15 MHz)^', consistent with each other within fit errors. Er- 
rors at low frequency (high SNR) are primarily due to systematic 
uncertainty in the antenna gains, and to thermal noise statistics at 
higher frequencies. 



originates in the earth's atmosphere and which involves elec- 
trical current accelerating transverse to the geomagnetic field. 
Such observations are in every way consistent with predic- 
tions of geosynchrotron emission from cosmic-ray air show- 
ers. In addition, the inherent spectral and time-domain simi- 
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FIG. 4: Plane of polarization of UHECR events compared to the an- 
gle of the magnetic field local to the event, with the red line indicating 
the expectation for the Lorentz force. The reflected events are cor- 
rected for their surface Fresnel coefficients, and angles are measured 
from the horizontal. 



larity of our radio pulses, as well as their robust correlation to 
geomagnetic parameters, suggests that ANITA's observations, 
which are at much greater distance and higher frequency than 
prior and current air-shower geosynchrotron observations, are 
less susceptible to near-field fluctuations of radio strength and 
plane of polarization. Such issues have been problematic in 
this field throughout most of its history. 

Our data represent the first broadband measurements of 
geosynchrotron emission in the UHF frequency range. The 
average observed radio-frequency spectral flux density of the 
above- and below-horizon events, shown in Fig.[3](Bottom) is 
consistent with an exponential decrease with frequency. The 
lack of any statistically significant difference in the spectra 
for the direct and reflected events indicates that ice rough- 
ness is unimportant for the average surface reflection. To es- 
timate the electric field amplitude at the source of these emis- 
sions, we model the surface reflection using standard physical- 
optics treatments developed for synthetic-aperture radar anal- 
ysis. Such models use self-affine fractal surface parame- 
ters [23] and Huygens-Fresnel integration over the specular 
reflection region to estimate both amplitude loss and phase 
distortion from residual slopes or roughness. In our case, 
we used digital-elevation models from Radarsat |'24'] to esti- 
mate surface parameters for each of the event reflection points, 
known to a few km precision. In most cases the surface pa- 
rameters are found to be smooth, yielding only modest effects 
on the reflection amplitude; in a minority of the events, sur- 
face parameters were estimated to be rougher, but still within 
the quarter-wave-rms Rayleigh criterion for coherent reflec- 



tion ||25|] . Fresnel reflection coefficients were determined us- 
ing a mean near-surface index of refraction of n == 1.33, typi- 
cal of Antarctic firn. 

To estimate the primary energy for the observed events, we 
used two independent approaches that determine the ampli- 
tude of the radio emission and the mean angular offset of 
the observed events. One approach is based on current air- 
shower geosynchrotron radio emission simulations developed 
for surface arrays jlSi Il9l 12011 . and the second approach is 
based on data-driven maximum-likelihood modeling in which 
a small set of parameters of a semi-empirical model were it- 
eratively fit to the observed characteristics and total number 
of the events, given the known UHECR energy spectrum Ij]]. 
The former method had the advantage of extensive work done 
to develop full-scale air shower Monte Carlo simulations for 
such radio emission; however, the simulations are not directly 
relevant to the very different geometry and higher frequency 
range of ANlTA's observations, which are in the far-field com- 
pared to most ground array observations, and which also in- 
volve showers at much larger zenith angles than ground ar- 
rays usually observe. The latter data-driven approach used 
physically-motivated parameterizations to capture the radio 
emission characteristics. The resulting constraints imposed by 
the data, including amplitude, phase, and frequency-spectral 
content for all 16 events, were found to be effective in fit- 
ting both the primary energy and observed angular distribu- 
tion of the events. Using the ground-based geosynchrotron 
models we found no self-consistent solution for the event en- 
ergy and mean angular offset, and we conclude that ANlTA's 
observations are in tension with the current ground-based sim- 
ulations. Our data-driven approach converged on a solution 
which gave estimated event energies as shown in Fig.|5]along 
with an overlay of the histogram of the energy distribution 
of simulated events seen in reflection. The implications of 
the data-driven solution are that the RF signals from these 
highly inclined, distant showers are significantly stronger than 
predicted by current geosynchrotron models. The mean en- 
ergy of the ensemble of reflected events is estimated to be 
1.5 ±0.4{stat)tl'i,isys) x lO'^ eV, approaching flie thresh- 
old of the Greisen-Zatsepin-Kuzmin (GZK) cutoff 1261 uTT . 
which marks the beginning of the absorption edge of UHE- 
CRs against the cosmic microwave background radiation. The 
large asymmetry in the systematic uncertainty is due to the 
uncertainty in the angular offset, which tends to strongly bias 
toward underestimating the event energy in our models. For 
the direct events, the mean energy is lower due to stronger di- 
rect signals, but the acceptance - limited to a narrow angular 
band around the horizon - is also much lower 

Based on our data-driven semi-empirical approach, we es- 
timate a mean angle of observation relative to the true shower 
axis of (1 .5 ± 0.5)°. This angular precision is comparable to 
that of ground-based cosmic-ray observatories, and adequate 
to allow us to map these events back to the sky. The final 
error circle is 2° in diameter after convolving with angular 
reconstruction precision and the modest tilts of each event 
locale, determined from Radarsat images at 200 m resolu- 
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FIG. 5: Top: energies of detected UHECR events, with reflected 
events in red, direct (above-horizon) events in black, and the sim- 
ulated event sample (reflected events only) shown in blue. Bottom: 
Map in celestial (a, 8) coordinates of the ANITA events (circles) with 
2.0 degree radii, and nearby AGN (grey diamonds) from the Veron- 
Cetty catalog. The approximate energy for each event is color coded 
by the logjo of the estimated event energy. ANlTA's exposure is ap- 
proximately uniform across the band 5° > 5 > —30°. 



tion i24ll . The resulting map is shown in Fig. |5] Our event 
positions are uncorrected to the sky positions of the Auger 
Observatory UHECR events, and the ensemble is also uncor- 
rected to AGN in the nearby universe. This is expected for 
events in this energy range where intergalactic magnetic de- 
flection is significant. While our sample of UHE events is 
significantly smaller than the current totals for the Auger Ob- 
servatory IjI, according to our models the acceptance of this 
method of UHE detection continues to increase at high en- 
ergies, even beyond 10^" eV, whereas the acceptance of all 
ground-based UHECR observatories saturate well before this. 
Estimates from our simulations indicate that, after optimiza- 
tion for UHECR observation, a new 30 day flight of ANITA 
could detect a total of several hundred geosynchrotron events, 
with 60-80 above lO'^ eV, and -- 10 above the nominal GZK 
cutoff energy. We conclude that a balloon-borne observatory 
is viable at the highest cosmic-ray energies, and if the fidelity 
of models of the geosynchrotron process continues to improve 
at the rate it has in recent years, such an approach will be able 
to further elucidate possible correlations in cosmic -ray origin 
directions as well as the shape of the endpoint of the UHECR 
energy spectrum. 

We are grateful to NASA, the US National Science Foun- 
dation, the US Dept. of Energy and the Columbia Scientific 
Balloon Facility for their generous support of these efforts. 
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